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Geodesic flows around charged black holes in two 
dimensions 

RAVI SHANKAR KUNIYAL^ • 
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Abstract We study the evolution of timelike geodesics 
for two dimensional black hole spacetimes arising in 
string theory and general theory of relativity by solv¬ 
ing the Raychaudhuri equation for expansion scalar as 
an initial value problem. The possibility of geodesic 
focusing/defocusing is then examined accordingly with 
different settings of black hole parameters. In view of 
the geodesic focusing/defocusing, the critical value of 
expansion scalar is also calculated in each case. The 
effect of the charge and the cosmological constant on 
the evolution of expansion scalar for timelike geodesics 
is discussed in detail. 

Keywords Geodesics; expansion scalar; Reissener- 
Nordstrdm and stringy black holes; geodesic focusing. 


1 Introduction 


The classical black hole (BH) solutions in general rela¬ 
tivity (GR) in four dimensions (4D) are endowed with 
new features, if the theory is modified by low energy 


string corrections in the action (Garfinkle et al. 

1991; 

Campbell et al. 

1990l IShapere et al. 

1991 1. Action 


for a two dimensional (2D) gravity can be obtained if 
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one introduces a scalar field (whose logarithm is called 
dilaton in the notions of string theory). Non-minimal 
coupling of the dilaton and axion fields with gravity and 
other fi elds is one of the key aspects o f these modified 
actions( Bekenstein 1972 : Ghase IQTOh . The dynamics 


of gravity for a generic higher dimensional theory may 
therefore be treated as an effective theory at lower di¬ 
mensions by the introduction of dilaton fields. 

As it is well known that if mass of the BH is large com¬ 
pared to the Planck mass then all vacuum solutions 
of Einstein gravity are basically the approximate so¬ 
lutions of low-energy string gravity. The charged BH 
solntions emerging in string theory in which the dila¬ 
ton field conples to the Yang-Mills field are fnndamen- 
tally different from the classical Reissner-Nordstrdm 


BHs (jMignemi and Stewart 1119931) in GR. 


dPoisson 2004 

: Wald 

1984: Hartle 200.8: 

19851: Chandrasekhar 

00 

.loshi h997h. t 


step may be to study the geodesic motion for given 
test particles in a given spacetime. There are vari¬ 
ous interesting aspects related to the geodesic motion 
around a give n 2D BH spacet ime, including the energy 

the varying constants 
the kinematics of time- 


distribnt ion (IVagenas 


200 . 


problem! E.G. Vagenas l200 


2009l Dasgupta et al. 2008: Dasgnpta et al. 

2009: 

Ghosh et al. I 2010l: IPernando I 2012: Pugliese 

2011) 

and geodesic deviation! Kolev et al. 20031 Unival et al. ^ 

2014: R.Unival et al. 2014). Farther the evolntion 

of the kinematical quantities (i.e. isotropic e 
sion, shear and rotation) along the geodesic flow 
governed bv the Ravchandhuri eanations !Kar 

xpan- 

are 

2008: 

Kar and sengupta IboOVllDadhich 1 2005: Ehlers 

20071). 


These equations result from relating the evolution of 
deviation between two neighbouring geodesics to the 
curvature of spacetime. 

The structure of the paper is as follows. A brief in- 
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troduction to the spacetimes considered is presented 
in the next section. In section III, the kinematics of 
deformation is discussed, while the exact solutions for 
geodesics and geodesic deviation in this 2D stringy BH 
backg round have already been studied in our previous 
work (|Unival et al. II20I4I1 . In 2D case, the Raychaud- 
huri equation is just the evolution equation for expan¬ 
sion scalar while the shear and rotation vanishes man¬ 
ifestly. The evolution equations are analytically solved 
and critical values for the expansion scalar is calculated 
for geodesic focusing/defocusing in each case. In sec¬ 
tion IV, the role of initial conditions on geodesic focus¬ 
ing/defocusing is analysed and compared with another 
charged 2D BH spacetime arising in GR. The results 
are summarized in section V along with the conclusions. 


2 The Model and Spacetimes 

The general spacetime for spherically reduced 2D BHs 
is given as, 

ds^ =—X{r) dt^ + X{r)~^ dr^, (1) 

for specific choice of X(r) emerging in GR and other 
alternative theories of gravity viz. string theory. To 
understand the dynamics of a test particle in the 
spacetime given in eq. the general structure of the 
geodesic equations corresponding to the above line ele¬ 
ment is given by, 

i + X (r) X(r)~^ r i = 0, 

X\r) - 1 x\r)X-^{r) P = 0, (2) 

where prime (/) and dot (•) represent the differentia¬ 
tion with respect to r and affine parameter (r) respec¬ 
tively. The general form of the first integral of the 
geodesic equations (using the constraint equation as 
= —1 for timelike geodesics) for above case 
can be obtained as, 


i = E X-^{r), 

r = ±^E^-X{r), (3) 

where ill is a constant. We are considering the following 
BH spacetimes for our investigation. 


Case I: The line element corresponding to the 2D 
charged stringy BH, arisi ng from heterotic string theory 
( McGuigan et al. 11199211 in the “Schwarzschild” gauge 
is given as. 


+ ql 


where, 0 < t < -l-oo, r_|_ < r < -|-oo, r+ is the event 
horizon of the BH. The parameters m, and A are the 
mass, electric charge and cosmological constants of the 
BH respectively and horizons are placed at, 


’'± = 2A*" 


-± - ql 


( 5 ) 


The asymptotic flat region is r = -|-oo (since q^ and 
A are positive) and the curvature singularity is at 
r = —oo. 

First integrals of geodesic equations given in eq.© re¬ 
duce to the following form for the above spacetime spec¬ 
ified by the eq. 

E 

1 — 2 m 6“^'’'*' -I- ( 7 g ’ 

r = ■sjE'^ — 1 + 2 — q'^ . ( 6 ) 


Case II: The met r ic (13 represents usual charged BH 
( McGuigan et al. 1199211 in 2D, i.e. the two dimen¬ 


sional analogue of 4D Reissner-Nordstrdm BH space- 
time with, 


V(r) = l-^ + ^, (7) 

where the asymptotic flat region is at r = oo, and the 
curvature singularity is at r = 0. Here m and q^,^ are 
the mass and electric charge of the BH. The positions 
of horizons are given as, 

2 ■ ( 8 ) 

First integrals of the geodesic equations ©, now reduce 
to the following form for spacetime given in eq. 0, 

E 

\ _ 2 m N ’ 

r r'^ 

r=±JE^-{l- — + ^). (9) 

There are two singularities present for both of the 
spacetimes given by eq. o and eq. 0. For extremal 
case in which charge parameter is equal to the mass, 
both singularities are merged, hence only one event 
horizon exists. It is worth noticing that the case where 
charge dominates over mass represents a naked singu¬ 
larity (i.e. a singularity without event horizon). 


X(r) = 1 — 2 TO e 


( 4 ) 
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3 Kinematics of Deformation 


3.1 The Raychaudhuri equation for expansion scalar 


Now the first approach to have an expression for the ex¬ 
pansion scalar (6) of timelike geodesics is to calculate 
0 as a function of r from the expression of the corre¬ 
sponding first integral equation u® = (t, f) and using 
6 = ViU^. The expression of 9 for spacetime given by 
eq.® using eqs. reads, 

^ _ X'{r) _ -2 TO A e-2^’'-f 2 

^ ^ ± ^— 1 + 2 m e~‘^^ ^ ~ Qs ^ 

(10) 

This expression for 0 diverges to infinity at a turning 
point of the geodesic motion (i.e. r = 0). But whether 
geodesics will focus or defocus depends on the implicit 
relation between different parameters given in numer¬ 
ator of the expression as well as on the sign of radial 
velocity, f. For r —> O’*" focusing of geodesics will oc¬ 
cur if TO A > while geodesics will defocus if 

TO A < q^. Similarly, for f —> 0“ one will have 
geodesic focusing with to A and defocusing 

with TO A > g^ 

While the expression of 9 for spacetime given by eq. © 
using eqs.® is given as, 

-TO r + ql^ 

±r2 r-2 - r2 + 2 TO r - g2^ ■ 

( 11 ) 

This expression for 9 also diverges to infinity at f = 0. 
Here again condition for geodesic focusing or defocusing 
is given by the implicit relation between different pa¬ 
rameters given in numerator of the expression as well 
as on the sign of f. For f —0+ focusing of geodesics 
will occur if TO r > g^^, while geodesics will defocus 
if TO r < g^„- Similarly, for r —>■ 0“ one will have 
geodesic focusing with m r < g^^ and defocusing with 
m r > g2 ^. An explicit dependence of 0 on r can be 
obtained by substituting r(r) from the solution of cor¬ 
responding geodesic equations for r given in eq. ® and 
eq.®. In the expressions of 9 given in eq. (fTUl) and 
eq. dlB, there is no way to specify initial conditions on 
9. Hence the effect of initial conditions on the evolu¬ 
tion of geodesic congruences cannot be discussed using 
these expressions. So in the next step to analyse the be¬ 
haviour of geodesics congruences on imposing the initial 
conditions we will consider the Raychaudhuri equation 
for expansion scalar along with the geodesic equation 
derived for the above mentioned spacetimes. 


=- 


X\r) 

2 f 


The evolution equations for timelike geodesics for 2D 
BH spacetimes consist of the Raychaudhuri equation 
for expansion scalar (as shear and rotation both van¬ 
ish for 2D matrices ( Dasguota et al. 20091) 1 and the 
geodesic equation derived for a particular metric. The 
Raychaudhuri equation for expansion scalar fo r gen¬ 


eral 2D BH spacetime (IDasgupta et al. 
in eq. o is given by, 


200911 given 


9 + 9^ = = -RimEu^ 


( 12 ) 


The general form of ea. dT^ can also be written as, 

9 + 0 ^-^ = 0 (13) 


where R is the Ricci Scalar (see Appendix I). The gen¬ 
eral Raychaudhuri equation for expansion scalar given 
in eq. m reduces to the following forms for the space- 
times used for the present study: 

the general Raychudhuri equation given in eg. (1131) for 
the expansion scalar in the background of spacetime 
given by eq. o for E^ = 1 reads as. 


9 + 9^- 


8 TO^ A^ 


■-b 


32 A^ TO^ g^ 


4r2 A2 TO2 + g2 (4 t 2 A2 TO2+g2)2 


0 (14) 


this equation will now be solved analytically (in next 
section) for the expansion scalar in view of different 
charge to mass ratios as, g^ >> m, q^ = to and q^ << 
TO (where to and q^ correspond to mass and charge of 
the stringy BH). 

The general Raychudhuri equation given in ea. (fT51) for 
the expansion scalar in the background of spacetime 
given by eq. o for E^ = 1 reads as. 


9+9^-- 


16 TO^ 


--b- 


48 gL 




= 0 


(15) 


this equation will now be solved numerically (in 
next section) for the expansion scalar with conditions 
Urn >> m, qRN = ™ and q^^^ « m (where m and 
g^jjy correspond to mass and charge of the Reissner- 
Nordstrdm BH) and p = (3 r -b yjq^j^ -b 9 t2 to"^)^/^. 


3.1.1 Evolution of expansion scalar for different cases 


In the present section corresponding Raychaudhuri 
equation m for expansion scalar 9 is solved analyt¬ 
ically, as an initial value problem with different condi¬ 
tions on BH parameters. 
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(a) Charge dominance over Mass i.e. Naked 
Singularity {q^ >> to); 

The solution of ea. (fT4)l is given by, 

. __ g Cl - 32 TO^ qj T _ 

(/3 Cl — 12 TO^ + 9 s)(4 A^ to^ + qI) 

(16) 

where a = 16 A"' + 24 A^ q'^ “ 3 

P = 4: X'^ m'^ — 3 q'^ T. 

and Cl is the integration constant written in terms 
of the initial conditions as, 

_ A 00 - 32 TO^ A^ To 

^ i? 00 — 16 Tq — 24 to^ A^ Tq + 3 g| 

(17) 

where, A = 48 to^ A^ Tq + 8 A^ Tq — q'^, 

B = 16 TO^ A® Tq — 8 TO^ g^ A^ Tq — 3 g^ ro. 

(b) Extremal Black Hole Case (g^ = to); 

The solution of ea. ifHl) is given by, 

128 C2 + 64 C2 - 12 r2 + 72 Cz r + 3 
~ (16 C2 t4 + 24 C 2 r2 + 3 r - 3 C2)(4 r2 + 1) 

(18) 

where C2 is the integration constant given in terms 
of the initial conditions as, 

_ 3 (4 00 A^ to 2 T,f - 4 A2 to^ + 00 g^ tq - qj) 
^ F 00 - 128 773 - 64 772 g| r3 - 72 77 g^ tq 

(19) 


where, 

F = 64A®TO®ro + 112A'^TO^g^ro + 12A2TO^ggro — 3g®, 

77 = A 2 to^. 

(c) Mass dominance over Charge {q^ << to); 
The solution of ea. flS is now given as, 


0(t) 


4{2mXTtan ^(^^^) + 2 C 3 ATOT + ggAm) 




( 20 ) 


where, k = (4A^TO2r2 + q^)tan 1 ( 21 ^ 12 ) ami 

^ = 4 C 3 A^TO^T^ + 2ATOggT + C 3 gg. C 3 is the integra¬ 
tion constant given in terms of the initial conditions 
as. 


C 3 


H6q-\- Jtq — Am A gg 
4 00 A2 m? Tq — 8to2 A2 to -I- 0o g^ 


( 21 ) 


where, H = (^AX^ m'^ Tq + q1) tan ^(-^-^), 

7s 

J = 2 00 A TO g — 8 to^ A^ ). 

7s 

where ro and 0 o are the initial values of affine parameter 
and expansion scalar respectively. 


4 Analysis of geodesic focusing and defocusing 

The occurrence of finite time singularity can be ob¬ 
tained from the exact solution of expansion scalar 
by solving the corresponding Raychaudhuri equation. 
From ea. dTB)) (i.e. the solution of Raychaudhuri equa¬ 
tion for gg >> to), the critical value of expansion scalar 
can be expressed in terms of initial values of various BH 
parameters as, 

16 Cl A"^ TO^ + 2AK to 2 g2 a 2 To - 3 Ci g^ 

=- LN - 

( 22 ) 


where. 


F=(CiTo-^), 


L = 4A^ TO^ '^0 + 7s, 


A = 4 Cl TO^ A^ Tq — 12 A^ TO^ Tg 


3 Cl gg To -bgg. 

(23) 


It can be observed from ea. d^^ that 0 ^ 00 as A —>■ 0, 
hence Ci can be fixed from this condition as. 


12 Tq A 2 to^ — g^ 

4 r)) A2 777,2 - 3 To g| ‘ 


(24) 


For a particular set of parameters used in eg. (1241) (with 
gg = 100, TO = 1, A = 1, To = 20), Cl is calculated 
and used in ea. (l2^ to obtain the critical value of ex¬ 
pansion scalar as, 0(j(3OO) = —0.02122. For 0o < 0g, 
geodesic focusing occurs (see fig.(II|)). A similar analy¬ 
sis for TO = gg (with gg = 1, to = 1, A = 1), results in 
0fj(O.5) = —1.65517, while for to >> q^ (with q^ = 1, 
TO = 100, A = 1) 0o(l) = 2. For to >> q^ geodesic 
defocusing is observed if 9q > 9q, while geodesic focus¬ 
ing does not occur otherwise. Hence it is clear that 
evolution of geodesics depends on initial conditions of 
expansion scalar critically. 

In fig.(P geodesic focusing is observed for gg >> to, 
hRN >> ITT- and q^ = m, q^^ = to cases, while this effect 
is absent if to >> gg, m >> Fig.([T](f)) shows that 

in the absence of charge g^^^, dotted curve represents 
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the behaviour of an analogue to spherically reduced 2D 
Schwarzschild BH spacetime. Due to stringy correc¬ 
tions made in the action, corresponding charged stringy 
BH spacetime does not reduced to Schwarzschild BH 
spacetime in the absence of the charge q^. Fig.Q rep¬ 
resents quantitative effect of variation of BH param¬ 
eters such as charge and mass of stringy BH on the 
evolution of expansion scalar. Variation of parameters 
does not effect the evolution of expansion scalar quali¬ 
tatively, but only the critical value of expansion scalar 
for focusing(/defocusing) is changed correspondingly. 
Fig.([3|) represents the effect of variation of the cosmo¬ 
logical constant A on the evolution of expansion scalar, 
it is observed for » m and q^ = m that geodesics 
focus at smaller values of affine parameter as A is in¬ 
creased. Variation of A does not effect the evolution if 


m»qg. 

The variation of time to singularity! Daseunta et al~l 
2 OO 9 II (Tg) versus initial value of expansion scalar (0o) 


for stringy BH is presented in Fig.®. Time to singu¬ 
larity reduces as mass dominates over charge, again this 
figure also justifies the absence of focusing for m >> q^. 








Fig. 1 (a) gg = = 100, m = 1, 6lg(300) = -0.02122; 

(>3) Is = = 2, m = 1 ; (c) = m = 1, 

0§(O.5) = -1.65517; (d) q^ = = 1, m = 100, 0S(1) = 2; 

(e) Qs = few = 1, "i = 2; (f) gg = 1, q^^ = 0 (this case 
corresponds to 2D Schwarzschild BH), m — 100; with A = 1, 
where solid and dotted curves represent 6{t) for stringy and 
Reissener Nordstrom BHs respectively. 


5 Discussion and Conclusion 


In the present work, we have studied the evolution of 
timelike geodesics in the background of charged 2D 
stringy and Reissener-Nordstrdm BH spacetimes. The 
conclusions are summarized as below; 

• The evolution of expansion scalar for a charged 
stringy BH and Reissener-Nordstrdm BH arising in 
GR in 2D is discussed in detail. The expression of 
expansion scalar and the conditions for occurrence 
of geodesic focusing/defocusing are examined using 
corresponding Raychaudhuri equations for expansion 
scalar with different settings of BH parameters. The 
effect of charge and cosmological constant on the evo¬ 
lution of geodesic congruences is also analysed. 

• Geodesic focusing is observed in the cases of naked 
singularity (i.e. when charge dominates over mass) 
and extremal BHs (i.e. when charge and mass are 
equal), while it is observed that focusing is ab¬ 
sent when mass dominates over charge for both the 
cases. The results obtained for stringy BH case are 
i n agreement w i th th e results of our previous work 


( TJnival et al. 2014[l in view of the solutions of 


geodesic deviation equation. 

By solving Raychaudhuri equation for expansion 
scalar as an initial value problem, it is clear that ini¬ 
tial values of expansion scalar plays a critical role 
in geodesic focusing/defocusing (wherever observed) 
in each case. It is observed from the evolution of 
expansion scalar that initially converging geodesics 
converge more rapidly in comparison to initially di¬ 
verging ones. The critical values of expansion scalar 
for stringy BH case are also calculated for q^ » m, 
q^ = m and q^ « m using different values for con¬ 
cerned parameters. It is observed that for q^ » m 
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(a) 





-\-2 -X-4 - X-6 


I-X-2-X-4-XHi 

(b) 



Fig. 3 Quantitative effect of variation of cosmological con¬ 
stant for stringy Bfl with,(a) = 100, m = 1, A = 1; (b) 

= m = 1, X = 1- (c) = 1, m = 100, A = 1. 


Fig. 2 Quantitative effect of variation of charge (in fig.(a)) 
and mass (in fig.(b)) for stringy BH; (a) m = 1, A = 1, ( 
figure at right side of fig. (a) with no number represents the 
zoom out section of fig.(a) near r = 0 ); (b) = 1, A = 1. 

and Qg = m the geodesic focusing occurs if 0o < 
while for m >> geodesic defocusing is observed if 
00 > 

• The pattern of geodesic focusing is found qualita¬ 
tively similar for the cases » m and q^ = m 
for different values of the cosmological constant while 
geodesic focusing is not observed in case of m >> q^. 
However the scale of geodesic focusing for the above- 
mentioned cases differs quantitatively. Further a 
slightly positive increment in cosmological constant 
value assists the scale of geodesic focusing such that 
focusing occurs more rapidly. 

• One can observe for stringy BH background that 
time to singularity (rg) decreases considerably for 
extremal case (i.e. q^ = m) in comparison to naked 
singularity (i.e. q^ >> m). It is also clear from 
the variation of Tg with initial value of the expan¬ 
sion scalar (0o) that geodesic focusing is absent for 
m » qg case. The comparative study of 2D charged 
stringy BH and Reissener-Nordstrdm BH shows that 
the qualitative behaviour of the evolution of the ex¬ 
pansion scalar is similar in both the cases while it 
differs quantitatively. 


In order to have more accurate visualization of the kine¬ 
matics of geodesic flows in the background of differ¬ 
ent charged BH spacetimes, we intend to investigate 
the Reissener-Nordstrom BH spacetime in 4D and its 
charged analogue in other alternative theories of grav¬ 
ity in future. The investigations made in this work can 
play an important role in addressing the issues related 
with the formation of accretion disk around BHs and is 
a matter of separate discussion. It would also be worth 
studying the nature of null congruences in a similar 
fashion along with the notions of gravitational lensing. 
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Fig. 4 Variation of time to singularity (rg) with initial 
expansion scalar (So) for stringy BH with, (a) = 100, 

m = 1, A = 1; (b) Qg = m = 1, A = 1; (c) Qg = 1, m = 100, 
A = 1. 


Appendix I 

Non-zero components of the Christoffel symbols and 
Ricci scalar: 


Case I: 

p 2Ame“^'^’’ — 

Fi^o = -2A(-1 -f 

^ 2A(me“^'^’’ — ggC”"'^’’) 

11 

R = 8A2(me-2^’' - 2qle-‘^^^) 


Case II: 


— — 
^ 01 — 


mr — q 


r (-r^ -I- 2mr - q\^) 


Too = + )(TOr-g2 ) 


i 11 — 


mr — q 


r (-r^ -I- 2mr - q\^) 
i?= 4(2mr-3g2 ) 
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